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ABSTRACT 

The model presented in this thesis is a computer, simulation model 
of ASW interactions between a formation of high value group ships, 
protected by some screening ships, and some penetrating submarines. 
The model is designed for use as an aid in improving the ability of a 
proposed screening tactic in the detection of a penetrating submarine. 

A systematic procedure to improve a screen's effectiveness against a 
known submarine threat is demonstrated,, and an example problem is 
worked using this procedure. 
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I. INTRODUCTION 



In the present era of a continuing and growing Sino-Soviet sub- 
marine threat, the problem, of effective anti-submarine warfare (ASW) 
to protect a high value group of shipping is of considerable interest to 
the Naval Officer. In recent years a great deal of effort has been 
expended in an attempt to develop an optimal solution to the ASW prob- 
lem; but, due primarily to the complexities of the environment and the 
resulting inability to quantitize parameters,. no> such optimal solution 
has been achieved. It is not the purpose of this paper to provide 
another solution technique for the ASW problem, but rather to provide 
a computer simulation model of an ASW action and to suggest one pos- 
sible use for the model. 

The model presented in this paper is designed for use as an aid in 
the design of ASW screens by providing measures of effectiveness of a 
screen in the detection of penetrating submarines, the number of attacks 
made by penetrating submarines, and the number of successful attacks 
made by a penetrating submarine. Using the computer simulation 
model, the effectiveness of a simple ASW screen against penetrating 
submarines is maximized within the constraint limitations provided. 
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n. DESCRIPTION 



The model is a computer simulation of the ASW environment, 
utilizing an event store procedure. The model is written in FORTRAN 
IV for the IBM 360/67 of the U. S. Naval Postgraduate School. A copy 
of the program is contained in Appendix I. 

Before beginning the model description, it is desirable to define 
a few terms commonly used in the model. 

A replication is defined as one complete ASW action. The action 
begins with the initial courses, speeds, and locations set by the user 
for the ships and submarine and proceeds until termination of the 
action. 

A run is a collection of one or more replications, each of which 
has the same initial values of initial courses, speeds, and locations 
of ships and submarine. The variation among replications in any run 
are caused by the randomly generated occurances within each replica- 
tion. 

A high value group is a group of one or more ships to be protected 
by escort ships from submarine attack. 

A screen is defined as a barrier formed by a group of escort ships, 
distributed about the high value group, whose purpose is to protect the 
high value group from submarine attack. In both actuality and the 
model there are two general screen types, the random patrol area 
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screen and the fixed station screen. In the random patrol area screen, 
screening ships are assigned a circular patrol area with a specified 
radius and a center that is fixed at a relative range and bearing from 
the high value group. The screening ships randomly patrol this area 
at some speed greater than that of the high value group. In the fixed 
station screen, the screening ships are assigned a fixed station relative 
to the high value group and maintain this station throughout the 
replication. 

CPA is defined as the closest point of approach of a submarine to 
a screening ship. 

Detection is defined as occurring when a penetrating submarine 
reaches a point within the detection range of a screening ship. 

A formation is defined as the group of ships that make up both the 
high value group and its screen. The center of a formation is referred 
to as ZZ and may or may not be occupied by a ship. 

The m.odel is designed to simulate an ASW environment. The 
scenario consists of a formation, proceeding along a predetermined 
track and being attacked by submarines attempting to penetrate the 
screen in specific places. The entire ASW action takes place suffi- 
ciently far in the first quadrant of an X-Y grid to ensure positive 
locations at all times in the model. The ability of the individual 
screen ship to detect the penetrating submarines is determined by a 
detection range that is uniquely determined for each screen ship in 
each replication. If a random patrol area screen is used, the ability 
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of a screen ship to detect a submarine is also influenced by the location 
of the screen ships within their patrol area;’ i;..e..,, at the extremities of 
the area, the screening ship may be several thousand yards closer or 
farther away from a penetrating submarine than, ta the center of the 
ship's patrol area. To use the model, the user must supply certain 
data consisting primarily of the initial locations,, courses and speeds 
of all ships, and some cui*ve of detection range vs., probability of 
detection for the individual ship's sensors. The specific data are 
placed on control cards in a standard format and, a. new set of control 
cards is u£;ed in each run. The flow of simulation log is contained in 
Figure 1. 

There are five events in the model. These are: 

1. A course change by a screening ship. This event, called 
whenever a screening ship reaches the edge of its patrol area, gen- 
erates a new course for the screen ship to continue its random patrol- 
ling. 

2. A detection by a screening ship. In this event the location of 
the detecting screen ship, detected submarine, and ZZ are noted at 
the time of detection. 

3. A CPA by a screening ship. At CPA,, the time; CPA range; 
screen ship and submarine identitites; and locations of screen ship, 
submarine, and ZZ are stored in a vector for reproduction at the end 
of the run. 



8 




Figure 1: Simulation Logic Flow 
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4. An attack by the submarine on a high value group ship. The 
time and location of the submarine and attacked ship for each attack 
are noted in this event. 

5. Termination of the replication. This event locates all screen- 
ing ships, submarines, and high value group ships at replication ter- 
mination. The mode of termination, either by successful submarine 
penetration or by a sufficient number of screen ship detections on. the 
submarine, is indicated. 

The measures of effectiveness for a specified formation computed 
in the model are: 

1. The number of attacks on high value group ships made by 
penetrating submarines. 

2. The number of times the penetrating submarines are detected 
by screen ships. 

3. The number of successful penetrations made by the submarine. 

These measures may be used to measure the effectiveness of a 

screen. The expression "effectiveness of the screen" refers to the 
probability of the screen ships detecting the penetrating submarines 
and, at the same time, detecting the submarine sufficiently far from 
ZZ to prevent attacks on the high value group ships. Clearly, a screen 
may be designed that will maximize the probability of detection of the 
penetrating submarine; but, with the necessity of minimizing the 
number of attacks made on high value group ships as an added con- 
straint, the problem of designing an effective screen is complex. 
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With the measures of effectiveness computed in the model, the user may 
evaluate the relative effect of changes in parameters; e.g^, changing 
the radius of the random patrol area on the effectiveness of the screen. 

It should be noted that the model docs not consider the aspect of 
screen ship attacks on detected submarines. There are two reasons 
for this. The first reason is that the inclusion of screeu ship attacks 
on detected submarines is not considered germaine to the measurement 
of screen effectiveness as previously defined. The second reason is 
that the inclusion of screen ship attacks would greatly complicate the 
model, adding to both time and size requirements. 
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UI. LIMITATIONS OF THE MODEL 



The limitations of the model fall into the areas of dimension limita- 
tions and ship maneuvering limitations. A more complete description 
of these limitations is contained in the following paragraphs. 

A. DIMENSION LIMITATIONS 

The model is limited to certain maximxim numbers of ships and 
replications in each run to preclude using excessive computer storage. 
These limitations are; 



1. 


Screen ships 


10 


2. 


Submarines 


10 


3. 


High value group ships 


10 


4. 


Replications 


20 



B, MANEUVERING LIMITATIONS 

1. Screening ships are constrained at all times to remain within 
their assigned screening stations. In the event of a detection of a 
submarine, the screen ships do not attack but remain in their stations. 

2. Submarine maneuvering is confined to determination of a 
course to pass under ZZ. 

3. High value group ships do not maneuver but maintain course 
and speed throughout the replication. 
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IV. ASSUMPTIONS 



In the model certain assumptions have been made to reduce the 
complexity of the simulation. These assumptions are discussed in the 
following paragraphs with respect to screening ships, penetrating 
submarines, and high value group ships. These assumptions are con- 
sistent with the philosophy of model design and the measures of effec- 
tiveness derived for the model. 

A. SCREENING SHIP ASSUMPTIONS 

1. Detection ranges are assumed to be Boolien in nature; detection 
occurs with probability one when the penetrating submarine is within 
this range, and detection does not occur otherwise.. The detection 
ranges are generated uniquely for each ship in each replication. 

Detection ranges are generated at the beginning of each replication 
using the following procedure: The model user provides some curve 

of detection probability vs. detection range for each sensor that is 
associated with the screen ships. A number is randomly drawn from 
a uniform distribution on [O, 1]; this number is the detection probability 
for the ship in that replication. The range corresponding to this detec- 
tion probability for the sensor is the detection range of the ship for that 
replication. 

2. As indicated above, the model is designed making an assumption 
of a uniform distribution on [O, l] for detection probability. If the user 
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feels that detection probabilities are distributed according to some dis- 
tribution other than uniform, the model subroutine generating screen 
ship detection ranges must be rewritten by the user- to reflect this 
different distribution. 

3. Screening ships remain within their assigned stations at all 
times in the model. 

4. Screen ships maintain a constant speed throughout each run. 

5. Upon detection of a penetrating submarine,, the detecting ship 
does not attack the submarine. The detection of the submarine is noted, 
and the replication continues. 

6. Once a screening ship has gained contact on a penetrating 
submarine, the interactions of that screen ship and that submarine 
are terminated; i. e. , as the replication continues,, the screen ship 
does not redetect the submarine or calcvilate a new CPA for the sub- 
marine with each new random course change of the screen ship in its 
patrol area. 

B. SUBMARINE ASSUMPTIONS 

1. The penetrating submarine conducts no evasive maneuvering 
in the model. This assumption corresponds to actual submarine 
maneuvering practices. The submarine fire control solution for a 
penetration yields a course to pass under ZZ., Elaborate submarine 
maneuvering in a screen penetration is avoided,, if possible, as any 
maneuvering has a detrimental effect on the submarine fire control 
solution. 
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2. The penetrating submarine conducts no evasive maneuvering 
while penetrating the screen, nor does it attack- screen ships. The 
submarine does conduct attacks on high value group ships when arriving 
within a weapon's release range set b-y- the- liser- for each run. These 
attacks allow records of the effects of submarine penetrations to be 
maintained. 

C. HIGH VALUE GROUP ASSUIVIPTIONS 

1. The high value group ships do not maneuver in the model but 
maintain the initial course and speeds assigned by the user for each 
run. 

2. High value group ship number 1 is always located at ZZ. If 
the formation for the high -value group does not have a ship at ZZ, high 
value group ship number 1 is a dummy at ZZ, and the maximum high 
value group size is reduced to nine. 
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V. INPUTS 



Inputs to the model are by control cards,, one set for each run. 

The entries on the control cards are all right adjusted; i. e. , the entry 
values set as far to the right as possible in the column width specified. 
The input cards are ordered sequentially as listed below and follow the 
source deck directly. 



A. CARD NUMBER 1 



The entry values of control card 1 are integer values. 



Columns 



Definition 



1-4 Number of screen ships in run 



5-8 Number of submarines in run 



9 -12 Number of high value group ships in run 
13-16 Number of replications in run 
17-20 Number of detections of submarine by screen 
ships to terminate replication 



Range 

I^N^IO 

I^N^IO 

I^N^IO 

1 ^ N ^^20 



l^N^ll 



21-24 Print screen ship detection range for each 
replication Si average detection range for 



each ship for the run 

25-28 Print out CPA Si locations of each ship Si 



0, No; 1, Yes 



submarine in replication 



0, No; 1, Yes 



29-32 Print out locations of screen ship, sub- 

marine, Si ZZ at time of a penetrating sub- 



marine detection by a screen ship 



0, No; 1, Yes 
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Col\amns 



Definition 



Range 



33-36 

37-40 



41-44 

45-48 

49-52 

53-56 

57-60 

6 1-64 



Minimum submarine penetration speed 
Print out locations of screen ships, sub- 
marines, &: ZZ at termination of each, 
replication 

Print out event calendar for each replication 
Print out screening ship information for each 
replication 

Print out high value group ship information 
for each replication 

Print out submarine information for each 
replication 

Submarine weapon release range (lOO's of 
yards ) 

Print out locations of high value group ship &r 
submarine at weapon release for each 
replication 



O^N 



0, No; 1, Yes 
0, No; 1, Y es 

0, No; 1, Yes 

0, No; 1, Yes 

0, No; 1, Yes 



0, No; 1, Yes 



B. CARDS FOR SCREEN SHIPS 

Cards for screen ships follow sequentially the first control card. 
These entries and those of following cards must be non-negative and 
may be integer or floating point. For each screen ship, there is one 
card with the following format: 
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Columns 



Definition 



1 - 8 


Screen ship number 


9 -16 


Initial course (in radians) 


17-24 


Ship speed 


25-32 


Initial X-position 


33-40 


Initial Y-position 


41-48 


Sensor type 


57-64 


Patrol station radius 


65-72 


Initial X-position of center of random patrol area 


73-80 


Initial Y-position of center of random patrol area 



C. CARDS FOR PENETRATING SUBMARINES 

Cards for penetrating submarines follow directly after the cards 
for the screen ships. For each submarine, there is one card with 
the following format: 

Columns Definition 



1 - 8 


Submarine number 


17-24 


Speed ( - minimum penetration speed) 


25-32 


Initial X-position 


33-40 


Initial Y-position 



D. CARDS FOR HIGH VALUE GROUP SHIPS 

Cards for the high value group ships follow directly after the 
cards for the penetrating submarines. For each high value group 
ship, there is one card with the following format: 
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Columns 



Definition 



1 - 8 


High value group ship numbei* 


9 -16 


Course 



17-24 Speed 



25-32 


Initial X-position 


33-40 


Initial Y-position 
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VI. OUTPUTS 



The outputs of the model are of two types,, general statistical out- 
puts and specific format outputs. The generaT statistical outputs are 
always provided by the model, while the specific format outputs are 
provided to the user only if they are desired-. Ta have specific format 
outputs, the user codes control card 1 in the appropriate column with 
a 1 for each specific format output desired and a 0 for those specific 
format outputs not desired. EIxamples of general statistical and 
specific format outputs are included in Figures 2,. 3, and 4. 

A. GENET.AL' STATISTICAL OUTPUTS 

1. The number of attacks on high value group ships by sub- 
marines for each replication and the total for the run. 

2. The number of detections made on the penetrating submarine 
by scr eening ships during each replication. The upper limit of this 
value is set on control card one for each run. 

3. The number of successful penetrations to ZZ made by penetra- 
ting submarines in each replication and the total for the run. 

4. Average detection range of the penetrating submarine from ZZ 
for each replication; average detection range and its standard deviation 
of the penetrating submarine from ZZ for each run. 
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B. SPECIFIC FORMAT OUTPUTS 



1 . Individual screen shlp^s generated detection ranges for each 
replication and the average generated detection range for the run. 

2. CPA of each screen ship to a submarine in each replication. 

3. Locations of ZZ, screen ship making detection, and the 
detected submarine at time of detection in each, replication. 

4. Locations of ZZ, screen ships, and submarines at the ter- 
mination of each replication. 

5. Event calendar for each replication.. 

6. Screen ship information at termination of each replication. 
This information consists of: 

a. Screen ship number 

b. Course (in radians) 

c. Speed 

d. Ship X- position 

e. Ship Y_position 

f. Sensor type 

g. Generated detection range 

h. Patrol station radius 

i. X-position of center of random patrol. area 

j. Y-position of center of random patrol area 

k. Time of next random course change 

l. Detection status of ship (1 if ship made detection, 

0 otherwise) 

m. Time of last position update 
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7. High value group ship information at termination of each 



replication. This information consists of: 

a. High value group ship number 

b. Course (in radians) 

c. Speed 

d. X_position 

e. Y-position 

f. Time of last position update 

g. Time first submarine passes under ZZ 

8. Submarine information at termination of each replication. 
This information consists of: 

a. Submarine number 

b. Course (in radians) 

c. Speed 

d. X_position 

e. Y_position 

f. Time of last position update 

g. Number of times the submarine has been detected by 
different screen ships. The upper limit of this value is set on control 
card 1 for each run. 

9. Locations of submarine and target ship at weapon release in 
each replication. 
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1. Examples of General Statistical Outputs 

a. Output for each replication 

WEPLAUNCH DETECTIONS PENETRATION 

1 1 L 

b. Output for each run 

PENETRATIONS WEPLAUNCH AVDET STDDEV 

8 12 6057.6 2386.65 

AVERAGE RANGE OF SUBMARINE DETECTION FROM ZZ 
IN EACH REPLICATION 

7765.18 8809.04 0.0 4286.71 8633..7S 6465.39 

2. Examples of Specific Format Outputs 

a. Individual screen ship's generated detection ranges for each 
replication & the average generated detection range for the run 
SHIP AVRANGE REP 1 REP 2 REP 3 REP 4 REP 5 

1 3049.9 0.0 1097.8 0.0 3380.6 7427.6 

b. CPA of each screen ship to a submarine in a replication 

SHIP SUB TIME CPA SHIPXLOC SHIPYLOC SUB XLOC 
3.0 1.0 73. 7245. 44238. 68293. 37006. 

SUB YLOC ZZ X LOC ZZ Y LOC 
67857. 40000. 64319. 

c. Location of ZZ, screen ship & detected submarine at time of 
detection 

DETECTION TIME 64 BY DD 4 ON SUBMARINE 1 AT 
X=36733.00 Y=62409.00 

Figure 2: Examples of Outputs Available in the Model 
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d. Locations of ZZ, screen ships, & submarines at replication 
termination 



TIME 


SCREEN SHIP 


X-POSIT 


Y-POSIT 


86 


I 


42180. 0 


7Z546.0 


TIME 


MAIN BODY SHIP 


X-POSIT 


Y-POSIT 


86 


I 


40000.0 


68652. 0 


TIME 


SUBMARINE 


X-POSIT 


Y-POSIT 


86 


I 


39910. 0 


68724. 0 


Event 


calendar for each replication 




TIME 


EVENT NO. 


SHIP SUB 


MAIN BODY 


86 


9 


1 I 


1 



f. Screen ship information at termination of each replication 
SmP COURSE SPD X- COORD Y- COORD SONAR 
4 1.62 15. 35434. 72029. 0. 



DETRANGE 


PA TRAD 


PATX 


PATY 


CCTIIvlE DS 


4691. 3 


1800. 


35830. 


71074. 


90. 1. 



XYTIME 

86 . 

g. High value group ship information at the termination of each 
replication 

SHIP COURSE SPD X- COORD Y-COORD XYTIME 
1 0.0 10. 40000. 68652. 86.0 

EVT9TIME 

86 . 

Figure 3; Examples of Outputs Available in the Model 
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h. Submarine information at the termination of each replication 



SHIP 


COURSE 


SPD 


X- COORD 


Y- COORD 


XYTIME 


1. 


1.28 


7. 0 


39910. 


68724. 


86. 0 



TIMES DETECTED 

1 . 

i. Locations of submarine & target ship at weapons' release in 
each replication 



TIME 


SUB NO 


SUB X 


SUB Y 


TARGET NO 


TAR X 


74 


1 


37230. 


67927. 


1 


40000. 


TAR Y 


LAUNCH RANGE 








64654. 


4000. 









Figure 4: Examples of Outputs Available in the Model 
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Vn. MODEL USE 



As an example of the use of the ASW simulation model, the prob- 
lem of detection for a screen of four ships is: optimized against a 
penetrating submarine threat. The use of the term, optimization in 
the context of this section refers to a search: for' an. improved, though 
not necessarily optimal, solution. The screen is specified to be a 
random patrol area screen with the four screen-ships symmetrically 
distributed about the formation bas e course in patroL areas whose 
centers are equally spaced. This symmetric distribution of screen 
ships will cover a specific arc in front of the high value group. The 
screen ships must be located a distance from ZZ that is at least 
equal to maximum submarine weapons' launch range. The submarine 
threat is known, with the submarine specified to pass between screen 
stations. The size of the main body is unspecified. (See Figure 5) 

A. PROCEDUT.E 

The procedure used for optimization in this thesis is that of 
response surface methodology (RSM). This procedure focuses on 
the class of simulations for which a single output variable exists, 
and the input and output variables are quantitative and approximately 
continuous . 
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— ^ . — Submarine penetration tracks 

r Radius of patrol area 

+ Center of patrol area 

^ ZZ 



Figure 5: Four Ship Random Patrol Area Screen with 

Specified Submarine Penetration Tracks 
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If there exists a functional relationship between the quantitative 
output variable and the quantitative input variables, this relationship 
may be expressed as: 



y = g (xj, . . . ,xj^) 



where: y denotes the quantitative response variable 

fir 

X. denotes the input variable 
g denotes the functional relationship 
and the value of the response variable y lies' on some k + 1 dimen- 
sional response surface. 

Assuming this functional relationship g to be continuous and 
differentiable, we will assume that a Ta.ylor series expansion of g 
can be expressed as: 

w k k 



y = g(X) = + :^ a.X. + ^ aj^XjXj 



+ . . . 

In the above notation, y is a measure of screen effectiveness deter- 



mined for a screen by an average of the number of submarine detec- 
tions over all replications in the runs; and is the value of the i^^ 
input parameter for the screen. 

The problem using the simulation model is that we do not know 
the functional relationship g, yet we wish to maximize g(X) subject 
to some set constraints; i. e. , maximize y subject to maximum and 

minimum values for X.. The model is used to estimate the coefficients 

1 

of the above Taylor series expansion by sampling. That is, the model 
will be used to obtain values for y for specified values of X^; and the 
technique of least squares can be used to estimate the values of the 
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coefficients of the Taylor series expansion. Having an approxinaation 
for the functional form of g(X), it is possible,, from this functional 
form, to estimate the gradient of g and to use this estimated gradient 
to determine the direction to move for fartheir sampling in an effort 
to maximize g(X). 

If the initial samples used to determine the values of the estimators 
of the coefficients of the Taylor series expajision are taken from a 
rather restricted area, it is reasonable to make the assumption that 
g is linear throughout this area. That is, the functional form of g(X) 
may be written: 

k 

y = gW = ^ a.X. 

and that the gradient of g is: 

(a 2 1 • • • t * 

The sampling procedure is to observe at least k + 1 values of y at 

k + 1 different values of X and to estimate the values of (a,-,ai, . . . , a, ) 

— U ^ k 

by least squares. A large number of simulation runs with different 
values of the X^'s may be made to determine the values of the a^'s by 
least squares. By considering only the largest and smallest values of 
the Xj's within the restricted area, there are 2 experiments that may 
be performed without redundancy to determine the value of the a^'s. 

The minimum number of experiments that will uniquely determine the 
values of the a^'s in k + 1 dimensional space is k + 1. Due to the high 
cost of computer time, the author did not consider it necessary to use 
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more than k + different values of the input variables to determine 

the value of the a.'s. 

1 

To ease the computational complexity of the solution, the user may- 
code each with +1 for its largest valu€^ and _L for its smallest value. 
B-y coding the input variables as +1, the assumed linear function 
relationship g is transformed to a new linear functional relationship 
of the form: 



y = b + 21 b.x. 
0 111 



where x. = _ ^ + 






A 



+ = maximum X- + minimum X- 



/\ = maximum X. _ minimum X. 

1 1 



and + £aj ^ + 

= aj 

By using k + 1 experiments, a unique solution for the estimated values 
of b^ through b^^ is obtained by the simultaneous solution of k + 1 
linear equations. 

The estimates of b^ through b^ obtained by sampling are then 

used as estimates of the gradient of the response surface. By using 

the gradient as a guide, one expects to proceed to an improved solution 

k 

V = b + w b.2 
o j 1 

where y denotes the predicted response along the path 
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(wb^, wb2, .... wb^) 

where w is a positive constant to indicate we are moving in the direction 
of increasing values of y. That is, we shall observe a new value of y 
at the point (wb^, wb^, . » . wb^^) and hope that the observed value of y 
will be close to the predicted value of y using y above. It should be 
noted that with increasing w, the accuracy of the predicted response 
may become unreliable unless the unknown functional relationship is, 
in fact, linear; and, thus, this path should be followed with caution. 

As long as observed values for y continue to increase along the path, 
one can continue. If the linear approximation for the unknown func- 
tional relationship has seemed satisfactory along the path, the appear- 
ance of stationary or decreasing values for y as the path is continued 
to be followed indicates that the linear approximation and gradient is 
calculated. The new path indicated by this gradient is followed until 
either the optimal solution is reached, or a constraint is reached; or 
a new linear approximation and gradient must be calculated. 

B. THE PROBLEM 

- For the problem being considered, the screen is specified to be a 
random patrol area screen with the ship located on bearing 300°, 340°, 
020°, and 060° relative to ZZ respectively. The screen center of 
each patrol area is to be no greater than 10, 000 yards from ZZ and 
no less than 4, 000 yards from ZZ. The submarine is to penetrate 
from 320°, 000°, and 040° relative to ZZ. The number of ships in the 
high value group is unspecified. (See Figure 5) 



31 



The model provides for each run the number of successful pene- 
trations made by the submarine. Therefore, using the number of 
unsuccessful penetrations as a measure of s creen effectiveness to 
optimize using RSM, one must first consider what variables in the 
model affect this measure of effectiveness. The variables that affect 
the output parameter of effectiveness y are; ship speed,, radius of 
random patrol area, the range of the center of the irandom patrol area 
from ZZ, ship detection range and the location of the ship within its 
random patrol area. The first three variables, ship speed, radius of 
patrol area, and the range of the center of the patrol area from ZZ, 
are variables over which the user has direct control in the model. 

Ship detection range and the location of the ship within its random 
patrol area are random within the model and, therefore, not directly 
controlled. 

For this problem, it was decided to limit the number of computer 
runs to eleven, with ten replications in each run, and to optimize the 
screen effectiveness using no more than those runs. Screen ship 
speed was fixed at 15 knots with the speed of the high value group at 
10 knots. Thus, there are two variables to consider in optimizing the 
screen: radius of the patrol area (X^) and the range of the center of 

the patrol area from ZZ (X^). 

For sampling, the initial area considered was bounded by the 
values of 1, 000 yards and 3, 000 yards for X^ and 8,. 000 yards and 
10, 000 yards for X 2 . The coded variables are; 
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■M 



X 

X = 1 - 2,000 

1 1,000 

X = ^2 - 9,000 

2 1,000 



The initial sample points and their associated r'esponses were: 



Sample point 




^2 


^2 


y 


1 


1, 000. 


10, 000 


-1 L 


0.533 


2 


1, 000 


8, 000 


-1 -I 


0. 567 


3 


3, 000 


8, 000 


-1 _L 


0. 667 



The solution of the three simultaneous linear equations yielded 
the following linear approximation; 

y = 0. 457 + 0. 017x^ _ 0. 075x2 

Since the largest observed value of the response was 0..667, it 
was decided to choose w such that the predicted response would be in 
the neighborhood of 0.700. The w used was 41, which corresponded 
to y = 0. 697 . We, therefore, sampled y at the point 
(41 (0. 017) , 41 (-0. 075)) 

The data for the next observation was as follows: 

Sample point w X. X x x y 

A ^ X ^ 

4 41 2,700 5925 w(0.017) w(-0.075) 0.667 0.697 

The value of the observed response for this run was 0. 667, a value 
less than predicted but still not less than the greatest previously 
observed value. It was decided to continue along this path to see if 
the solution would improve. Accordingly, a value for w of 49, 
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corresponding to a predicted response of 0.747, was chosen. We 
sampled y at the point 

(49 (0,017) , 49 (-0..075)).. 

The data for this next obs ervation was fallows: 

Sample point w y y 

5 49 2800 5325 49(0..0L7) 49(-0.075) 0.767 0.747 

By examination of the differences between pxedicted and response out- 
put for the previous runs, it appeared that llae linear approximation 
was not valid. Nevertheless, as long as the response continued to 
increase, the path indicated by our approximation of the gradient of g 
appeared to indicate the direction to an improved solution. Continuing 
along the path, we next selected w = 58, a value corresponding to 
y = 0. 800, and sampled at the point 

(58 (0. 017) , 58 (-0.017)). 

The data for this next observation was as follows: 

Sample point w X^ X^ x^^ x^ y 

6 ' 58 3000 4650 58(0.017) 58(-0.075) 0.667 0.800 

As the value of the response has -decreased, it was necessary to 
return to the point of maximum response and calculate a new linear 
relationship and gradient. Again, three sample points were used with 
the previous highest response point as one of these. These sample 
points and their associated responses were: 
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Sample point 




X 

2 




X 

2 


y 


7 


1800 


5325 


-I 


1 


0. 733 


8 


2800 


5325 


I 


1 


. 0. 733 


9 


2800 


4325 


I 


_1 


0. 800 



The solution of the simultaneous linear equations associated with 

these sample points yielded the following linear approximation: 

y = 0.7665 + Ox, - 0..0335x .. 

1 Z. 

It is of interest to note that the value o£ the estimator correspond- 
ing to the value of the radius of the patroXarea is zero. This indicates 
that the value of the radius is at its optimal value and can, therefore, 
no longer contribute to an improved solution. A value of x^^ =0 cor- 
responding to Xj = 2, 300 yards was, therefore, assigned in the final 
runs of the model. 

We decided to choose w = 30 such that the predicted response 

A 

would be y - 0. 800. We sampled y at the point 
(30 (0), 30 (-0.0335)). 

The data for this next observation was as follows: 

A 

Sample point X^ x^ • y y 

10 2300 4820 0 30(-0.0335) 0.800 0.800 

It was realized that the predicted response of sample point 10 was not 
an improvement over the best previous value; but the value of X^ was 
approaching the overall minimum constraint of 4, 000 yards, and a 
more cautiovis step wu.s taken. 
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Continuing along the path for the final run, a value of w = 74 
corresponding to y = 0. 850 was chosen. We sampled y at the point 
(74 (0), 74 (-0.0335)). 

The data for this final observation was as follaws:. 

A 

Sample point X X x x y y 

± ^ L ^ 

11 2300 4085 0 74(-0..0335) 0.833 0.850 

With the completion of the eleven runs of the model, the optimal 
screen for the given constraints was found to have values of Xj, and 
X^ of 2,300 yards and 4,085 yards respectively... The optimal response 
for these vs. lues of X and X was y = 0. 833; this means that the op- 

J. ^ 

timized screen can be expected to detect a penetrating submarine prior 
to the submarine attacking the high value group 83. 3% of the time. 

With additional runs, the response could possibLy be improved; but, 
with only the variable X^ available to vary and with its value so close 
to its lower limit of 4, 000 yards, the improvement would probably be 
marginal. A summary of the above steps is contained in Figure 6. 

In this sample example, a simple screen that had several con- 
straints was considered. By using RSM, the screen variables were 
optimized against the threat of submarine penetration. Using the 
same techniques, the model may be used to optimize more complex 
screens with additional variables. 
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Run No 



Y Comments 



^2 y 



1 


1, 000 


10, 000 


-1 


1 


0.533 


. Value for path of 
steepest ascent are: 


2 


1, 000 


8, 000 


-1 


1 


0.567 


(2000 + 17W, 9000-75w) 
iin. Xj^,X 2 space 


3 


3, 000 


8, 000 


1 


1 


0. 667 


y=0. 457+w(0. 005914) 


4 


2, 700 


5,925 


__ 


- 0.697 


0.667 


w = 41 




2, 800 


5, 325 


- - 


- 0.747 


0.767 


w = 49 point of 
maximum response 


6 


3, 000 


4, 650 




- 0.80 


0.667 


w = 58 


7 


1, 800 


5, 325 


-1 


1 


0.733 




8 


2, 800 


5, 325 


1 


1 


0.733 


y=0. 7665+w(0. 00112225 
b =0.0 in this solution 


9 


2, 800 


4, 325 


1 


1 


0. 80 




10 


2, 300 


4, 820 


_ _ 


- 0.8 


0.8 


w = 30 


11 


2, 300 


4, 085 


__ 


- 0. 85 


0. 833 


w = 74 



Figure 6: Summary of the Search for Improved Solution 
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VIII. CONCLUSIONS AND RECOMMENDATIONS 



The model presented in this thesis is a siraulatian. model of ASW 
interactions between a formation and penetrating submarines. A pro- 
cedure has been demonstrated in which the model was used to improve 
the effectiveness of a proposed screening tactic against a known sub- 
marine threat. 

An attempt to use the model to measure the effectiveness of 
various screen tactics was examined. The procedure used in this 
attempt was to determine the average detection range of penetrating 
submarines from ZZ for each screen and to compare these average 
ranges by the analysis of variance techniques. This procedure to rank 
screen effectiveness in order of average detection range from ZZ, 
with the screen having the greatest average detection range being the 
most effective, had intuitive appeal. Unfortunately, the standard 
deviations of the average detection ranges were so great, on the order 
of two to three times the differences between average detection ranges 
of the screens, that no statistical interferences could be made with 
this procedure. 

In an effort to use the same general procedure, an attempt was 
made to compare the effectiveness of proposed screen tactics by com- 
paring the average time until first detection of a penetrating submarine. 
This procedure also yielded no statistically significant results for 
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the same reason as the previous proposal. The results of these 
attempts are not included in this thesis. 

The model has not been successfully used to directly compare 
the effectiveness of various proposed screening: tactiics-; .but it may 
be used as shown in the preceding chapter to maximize the effective- 
ness of a screen, with certain constraints, againsta specific sub- 
marine threat. 
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APPENDIX A 



COMPUTER PROGRAM 
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33 RPLCSE=PI+THET A 

IF(BETA,GE,0.0.AMD. BETA. LE. (PI/ 2 . 0 ) )0n TO 80 
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EA3=ECAL( 1, 5) 

ECALd , 1 >=999999999 

RETURN 

END 
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SUBROUTINE FIND IS CALLED TO GENERATE TIME OF EITHER DETECTION OR 
CPA FOR THE SCREEN SHI^> ONLY IF THE SCREEN SHIP HAS NOT PREVIOUSLY 
MADE A DETECTION ON THE SUBMARINE. 
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STORE TIME OF NEXT COURSE CHANGE 
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EVENT CALENDAR AND DUMMY EVENT CALENDAR RESET FOR NEW REPLICATION 
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APPENDIX B: LIST OF PROGRAM VARIABLES 



This appendix contains the name and definitiorL of th.e symbols 
appearing in the computer program. This list contams only those 
variables using storage space in the program and does not contain 
those used only for temporary storage. 
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Name 



Definition 



CPA(I, J, K) 


CPA information for the ship onr. the submarine 



in a replication. The values hr this vector, are: 



CPA(I, J, 1) 


Screen ship number 


CPA(I, J, 2) 


Penetrating submarine number 


CPA(I, J, 3) 


Time of CPA 


CPA(I, J, 4) 


CPA range 


CPA(I, J, 5) 


Screen ship X location at CPA 


CPA(I, J, 6) 


Screen ship Y location at CPA. 


CPA(I, J, 7) 


Penetrating submarine X location. at. CPA 


CPA(I, J, 8) 


Penetrating submarine Y location at CPA 


CPA(I, J, 9) 


X location of ZZ at CPA 


CPA(I, J, 10) 


Y location of ZZ at CPA 


DECAL (I, J) 


Dummy event calendar, a copy of all events entered in 



the event calendar. The enti'y in DECAL contains: 



DECAL(I, 1) 


Event time 


DECAL(I, 2) 


Event number 


DECAL(I, 3) 


Number of screen ship in event 

/ 


DECAL(I, 4) 


Number of penetrating submarine in event 


DECAL(I, 5) 


Number of main body ship in event 


DESKIL 


An index of successful attacks by a penetrating sub- 




marine on a main body ship in a replication 
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Name 


Definition 


DSSV(I, J) 


Dummy submarine vector containing, initiail information 


DSSVd, 1) 


tVi 

Number of penetrating submarine 


DSSV(I, 3) 


Speed of submarine 


DSSV(I, 4) 


Initial X location of submarine 


DSSV(I, 5) 


Initial Y location of submarine 


DSTAT(I, J, K) 


tt 

Detection range from ZZ in the I replication for the 



submarine by the screen ship 



DSV(I, J) 


Dummy screen ship vector containing initial'- information 


DSV(I, 1) 


th 

Number of r screen ship 


DSV(I, 2) 


Initial screen ship course (in radians) 


DSV(I, 3) 


Screen ship speed 


DSV(I, 4) 


Initial screen ship X location 


DSV(I, 5) 


Initial screen ship Y location 


DSV(I, 6) 


Screen ship sonar type 


DSV(I, 8) 


Screen ship patrol area radius 


DSV(I, 9) 


Initial X location of center of assigned patrol area 


DSV(I, 10) 


Initial Y location of center of assigned patrol area 


DTARd, J) 


Dummy main body ship vector containing initial in- 



formation 



DTARd, 1) 


f "h 

Number of main body ship 


DTAR(I, 2) 


Initial main body ship course (in radians) 


DTAR(I, 3) 


Main body ship speed 


DTAR(I, 4) 


Initial X location of main body ship 


DTARd, 5) 


Initial Y location of main body ship 
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Name 


Definition 


EAl 


Number of a screen siiip 


EA2 


Number of a penetrating submarine 


EA3 


Nximber of a main body sbip 


ECAL(I, J) 


entry in event calendar 


ECAL(I, 1) 


Event time 


ECAL(I, 2) 


Event number 


ECAL(I, 3) 


Number of screen ship in event 


ECAL(I, 4) 


Number of penetrating submarine in. event 


ECAL(I, 5) 


Number of main body ship in event 


EN 


Event number 


ET 


Event time 


KT 


Index of total number of events entered’ in DECAL(I, J) 


MAXEV 


Maximum number of events allowed in ECALi(I, J) 


MBKIL 


Number of successful submarine penetrations in 
replication 


MBKILL 


Number of successful submarine penetrations in run 


MSTAT(I) 


Control vector for each run 


MSTAT(l) 


Number of screen ships in run 


MSTAT(2) 


Number of submarines in run 


lvlSTAT(3) 


Number of main body ships in run 


MSTAT(4) 


Number of replications in run 


MSTAT(5) 


Number of detections of submarine by screen ships to 
terminate run 
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Name 


Definition 


MSTAT(6) 


Print screen ship detection: range for:- each replication 



it average detectiaix range for- each-s-creen ship during 





run 


MSTAT(7) 


Print out CPA & locations of each sliip &: submarine in 



replication 



MSTAT(8) 


Print out locations of screen ship; submarine, and 



ZZ at time of detection of a penetrating submarine by 





a screen ship 


MSTAT(9) 


Minimum submarine penetration speed 


MSTAT(IO) 


Print out locations of screen ships,, submarines, & 



ZZ at termination of each replication 



MSTAT(ll) 


Print out event calendar for each. replication 


MSTAT(12) 


Print out ship vectors for each replication 


MSTAT(13) 


Print out main body vectors for each replication 


MSTAT(14) 


Print out submarine vectors for each replication 


MSTAT(15) 


Submarine weapon release range (lOO's of yards) 


MSTAT(16) 


Print out locations of main body ship submarine at 



weapon release for each replication 



MTARK 


Index of submarine attacks on main body ships 


PI 


3. 14159 


R 


Index used in subroutine FIND to isolate sections of 




the subroutine after the first call’ on the subroutine 
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Name 


Definition 


SONAR(I, J) 


The generated detection range of: the screen ship 



for the replication in a. run 



SONAV(I) 


The average o£ generated detection. ranges of the 



screen ship in a run 



SSV(I, J) 


til 

Information vector of the B snbmarine in a replication 


SSV(I, 1) 


Number of the submarine 


SSV(I, 2) 


Submarine course (in radians) 


SSV(I, 3) 


Submarine speed 


SSV(I, 4) 


Submarine X location 


SSV(I, 5) 


Submarine Y location 


SSV(I, 6) 


Time X-Y location fixed 


SSV(I, 7) 


Number of times detected by screen ships 


SUBKIL 


Number of times the penetrating submarines detected 


SV(I, J) 


til 

Information vector of the screen ship in a repli- 



cation 



sv(i, 1 ) 


til 

Number of the I screen ship 


SV(I, 2) 


Screen ship course (in radians) 


SV(I. 3) 


Screen ship speed 


SV(1,4) 


Screen ship X location 


SV(1,5) 


Screen ship Y location 


SV(I,6) 


Screen ship sonar type 


SV(I,7) 


Screen ship generated detection range 


SV(I, 8) 


Patrol area radius 
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Name 


Definition 


SV(I, 9) 


X location of center of assigned patmL area 


SV(I, 10) 


Y location of center of assigned patrol ar-ea. 


SV(I, 11) 


Time of next course change 


SV(I, 12) 


Number of submarines detected in replication 


SV(I, 13) 


Time X-Y location fixed 


TAR(I, J) 


til 

Information vector for the I main body ship 


TAR(I, 1) 


Number of the main body ship 


TAR(I, 2) 


Main body ship course (in radians) 


TAR(I, 3) 


Main body ship speed 


TAR(I, 4) 


Main body ship X location 


TAR(I, 5) 


Main body ship Y location 


TEMP(I) 


Used to aid in ordering event calendar as a temporary 




storage location for an event 
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